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The spectral weights (SW's) for adding and removing an electron of the Gutzwiller projected 
d-wave superconducting (SC) state of the t- J-type models are studied numerically on finite lattices. 
Restrict to the uniform system but treat exactly the strong correlation between electrons, we show 
that the product of weights is equal to the pairing amplitude squared, same as in the weakly coupled 
case. In addition, we derive a rigorous relation of SW with doping in the electron doped system 
and obtain particle-hole asymmetry of the conductance-proportional quantity within the SC gap 
energy and, also, the anti-correlation between gap sizes and peak heights observed in tunneling 
spectroscopy on high T c cuprates. 

PACS numbers: 



The emergence of the superconductivity as holes doped 
into the Mott insulating parent compounds is one of the 
intriguing phenomena of high T c cuprates 1 . It is usu- 
ally emphasized, however, below the transition temper- 
ature there is strong similarity of the superconducting 
(SC) state with that of low T c materials in the sense 
that the state is composed of Cooper pairs of electrons, 
though with unconventional rf-wave pairing symmetry 2 
and thin superfluid density 3 . Nevertheless, differences 
from the conventional SC state are unearthed clearly by 
high-resolution scanning tunneling microscopy and spec- 
troscopy (STM/STS) on different cuprates with easily- 
cleaved surfaces. Namely, despite of physical quanti- 
ties with nanometer scale inhomogeniety, the averaged 
(over some area within the scanned field-of-view) tun- 
neling conductance is of unexpected behavior: its am- 
plitude at negative sample bias- for removing electrons 
- is often larger than that at positive one- for adding 
electrons Most intriguingly, the gap size, which 

was usually inferred from photoemission experiments 8 , 
can be extracted directly from conductance peaks for the 
first time and is found to be larger as peaks become less 
pronounced^iSi Same behaviors have been observed pre- 
viously in underdoped cuprates using point-contact tun- 
neling set-up^. 

The theoretical attempt to understand the features 
seen in STM/STS hitherto is mostly about the effects 
of the inhomogeneous dopant induced SC order parame- 
ters to the local density of states (LDOS) of the d-wave 
BCS (d-BCS) stated, it largely neglects the strong cor- 
relation between electrons which should be essential for 
the case of underdoped cupratesii. Actually, there is not 
enough understanding about the effect of strong corre- 
lation for a homogeneous system to help us to address 
the complex issues of disorder as revealed by tunnel- 
ing experiments. Recently, Anderson proposed that the 
asymmetric tunneling conductance is closely related to 
the strong correlations inherent in the Gutzwiller pro- 
jected d-BCS or, simply, resonating-valence-bond wave 



function (RVB WF) 12 . However, treating the projection 
only approximately by the usual scheme of the renormal- 
ized mean- field theory (RMFT), there have been con- 
troversy whether the asymmetry is accounted for by the 
coherent quasi-hole (QH) and -particle (QP) excitations 
of the projected state or rather by the incoherent part 
dictated by the spectral sum ruleiii^i^. Also, the cor- 
relation between gap sizes and peak heights has not yet 
been examined clearly from the strong correlation point 
of viewiii. 

In this paper, we defer the issue of inhomogeneity to 
later work and examine exactly the effects of strong cor- 
relation by numerically investigating the spectral weights 
(SW's) of the d-wave RVB (d-RVB) state on finite square 
lattices. With SW, particularly, for removing an 
electron (defined in Eq.Q) below), calculated, we ob- 
tained several results: (i) d-wave pairing amplitude 
squared is equal to the products of SW's, as it is ex- 
actly for weakly-coupled case; (ii) inspired by the hole 
doped case we focus mostly here, a rigorous relation of 
SW for removing an electron is derived for the electron 
doped case; (iii) the difference between Z^ a calculated 
exactly and by using RMFT becomes significant at low, 
finite doping of holes; (iv) the particle-hole asymmetry 
of the sum of low-energy SW's over momenta within 
some energy window near the Fermi level becomes more 
prominent with reduced dopingii and, meanwhile, (v) 
the heights of SW peaks decrease as the gap sizes in- 
crease. 

Let us start by defining the SW for adding (and re- 
moving) one electron we calculate, i.e. 

7+( -) = i (N e M-n i <Lm i N e ) i 2 

kCT (JV e |JV e )<JV e +(-)l|JV e +(-)l)' {> 
where, for momentum k, 

| N e + 1) = P d cl a | N e ) (2) 



2 



for the QP excitation, and 

\N e -l) = P d d k _„ | N e - 2) c 
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for the QH one which is also proportional to PdCua- \ N c )q. 
Here | N e ) is related to the trial WF of the projected 
electron-paired ground state in a uniform system, 
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The variationally optimized | N e ) we focus on in this 
paper is the d-RVB stated. With N e the total number 
of electrons, coefficient a q = Vq/uq = (Eq — e q )/A q in 
which Vq and it q are SC coherent factors, e q = — (cosq x + 
cos q v ) — /i — t' v cosq x cos q y — t' v (cos 2q x + cos 2q y ), A q 



A (cos q x — cos q y ) , and _E q 



A q . The operator 



Pd projects out the doubly-occupied sites in the system 
with a finite number of doped holes present. In addition 
to A„ and t' v and ty are the other two variational pa- 
rameters associated with the long-range hoppings in the 
t-t'-t"-J model Hamiltonian, H = — Pdtij(c\ a Cj^ + 

h.c.)p d +jj: <i,j>(^i ■ S 3 — jUiUj), where hopping am- 
plitude tij=t, t' ', and t" for sites i and j being the nearest- 
, next-nearest-, and the third-nearest-neighbors, respec- 
tively, Si the spin operator at site i, < i,j > means that 
the interaction between spins occurs only for the nearest- 
neighboring sites. 

Applying identities for projection operator, 



[cka,Pd}Pd = 0; 



P d c kCT [4 v ,P d ] - P d {^- T J2e t( - k '~ k > R -n h - a }P d (5) 



with Ri a the position vector of the i-th spin a in the 
lattice of size N and n ia = c\ a Ci Gl we can relate the Z^ a 
exactly to the momentum distribution function (MDF) 
nicer as 



i 
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where x is the density of doped holes and rik = (-^e | 

<L<*° I N e)/(Ne I N e ^±ZL. 

As a digression to electron doped (ED) case, it is 
straightforward to show that, applying the hole-particle 
transformation to Eq. © 19 , SW of removing an electron 
in ED system satisfies rigorously ^ kcr =nko- — (1 — x)/2. 
This relation may be verified in experiment. 

Back to the hole doped case, although there is no exact 
relation like this for ZZ a , we notice that Z£ a and ZZ a 
satisfy a relation, 
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FIG. 1: The SC pairing amplitude for d-RVB state as a func- 
tion of doping determined by the products of SW's using 
Eqs.Q and @. System size here is 12 x 12. Different symbols 
represent results obtained for different values of (t',t")/t, as 
indicated 



which can be proved straightforwardly by combining 
Eqs.© and JSJ. The matrix element Pk, which repre- 
sents the off-diagonal long-range order in the pairing cor- 
relation, is related to the <i-wave SC pairing amplitude 
or order parameter A op by 



A p — 



2 

N 



| cos k x — cos k v 



(8) 



With both the SW's computed numerically, we plot in 
FigH the doping dependence of A op which indeed has 
the dome-like shape similar to the T c versus doping de- 
termined experimentally. Actually, the peak positions 
shown in FigQ are almost the same as what have been 
obtained previously by studying the d-wave long-range 
pair-pair correlationi&S. Furthermore, with more holes 
doped into the system, just like the reduction of long- 
range correlation between electron pairs is induced by 
the change of the anti-nodal Fermi surface geometrjii&, 
the SC order parameter is decreased due to ZZ a with k 
near (tv, 0). Hence Eq.Q provides another way to evalu- 
ate the strength of the pairing amplitude. 

For the BCS theory without projection, we know 



H-(-) 



l(v£) and Eq.Q) is also exactly satisfied. For 
the strongly correlated t- J-type models, even though the 
same relation is followed in RMFT— , it is still surprising 
to find out that this relation is correct in the RVB state 
with projection rigorously obeyed. 

On the other hand, reminiscent of what have been 
argued previously by analytic approach 13 , we recognize 
that the strong correlation effects becomes apparent only 
in ZZ at low doping. The effects due to strong correla- 
tion are examined by comparing the coherent SW's aver- 
aged over all momenta, i.e. Z~ ve = ZZ^/N, and the 
incoherent part defined by the relation 

n ave = n ave ~ Z ave (9) 

obtained by exact treatment of the projection and by us- 
ing RMFT. Here n ave = £ k n k(T /7V is the average MDF 
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FIG. 2: The doping dependence of SW's for removing an elec- 
tron and average MDF of d-RVB state obtained numerically 
for the 12 x 12 lattice and by RMFT. Numerical and expected 
n a ve are denoted by the circles and the solid line, respectively. 
The squares (triangles), connected by dashed (dotted) line as 
the guide for the eye, are for exact results of Z~ ve \n^ e , ex- 
tracted using Eq.@ with n ave calculated numerically]. All 
solid symbols are results for (£' , f")/£=(-0.3,0.2) and empty 
ones for (0,0). The dashed and dotted lines without data 
points represent results by RMFT. 



FIG. 3: g(E) for d-RVB state versus excitation energy E 
for 12 x 12 (solid circle) and 20 x 20 lattices (empty square 
and triangle, for different dopings). The associated excitation 
gap positions are marked by arrows (see text). Inset: Same 
plots for d-BCS state with gap value 0.6 and x = 0.125. The 
data shown here are for different sizes, denoted by the same 
symbols as in the main figure, with parameters fi/t = —0.269 
(12 x 12) and (i/t = -0.272 (20 x 20), respectively. All data 
are obtained for (t',t")/t=(-0.3,0.2). 



which should always be equal to the electron density of 
the system. 

The exact results for the 12 x 12 lattice and that 
by RMFT are shown in Fig [21 The coherent part of 
Zr by RMFT is gtv^ with renormalization factor g t = 
2x/(l + x). Completing the momentum sum for the co- 
herence factor, the average result is x(l — x)/{l +x) and, 
thus, <™ c e = [(l-x) 2 ]/2(l+x) u , plotted in FigH (dashed 
and dotted lines, respectively) in comparison with the 
exact ones. As is shown there, while the numerical n ave 
(solid circles) is indeed equal to the electron density, the 
exact incoherent SW for removing an electron is less than 
the RMFT result. The difference becomes more signifi- 
cant as hole doping level is reduced. Interestingly, this 
behavior is independent of the (t',t")/t values (repre- 
sented by solid and empty symbols in Fig[2) which corre- 
spond to very different doping dependence of the Fermi 
surface shape and also the DOS. By contrast, the aver- 
age values of Z£ calculated exactly (not shown) and by 
RMFT are identical due to Eq.©. 

To make a comparison with tunneling experiments, we 
then concentrate on the SW's as a function of the excited- 
state energy. By applying the model Hamiltonian to exci- 
tations | N e ±l), we calculate their excitation energies for 
each momentum and also the corresponding energy gap 
by fitting the excitation energy E^. To reduce the effect 
of finite size, we define the sum of Z^ a /N, over momen- 
tum k which has energy within E — AE/2 and E+AE/2, 
as g(E) [negative(positive) for removing(adding) an elec- 
tron] which could be viewed, approximately, as propor- 
tional to the conductance at low energy E. We plot g(E) 
in Fig|SJ up to about the energy where peaks appear 
for lattices of size 12 x 12 with AE/t =0.3, and also 



20 x 20 with energy interval 0.2 for various dopings 25 . 
To make sure our treatment is correct, we have also ap- 
plied the same analysis to the d-BCS state. As shown in 
the inset of FigEl the ideal BCS result is hardly distorted 
by the finite size. Note that, with the reasonable finite- 
size dependence, we obtain indeed the F-shape d-wave 
gap near zero energy. The width between peak posi- 
tions is also roughly equal to two times of the gap value 
deduced from the excitation energy. Looking at the re- 
sults closely, while g(E) may indeed be about the same 
at the opposite sides in the very vicinity of zero energy 
as suggested in Ref.22, g(E) for removing an electron is 
always larger than that for adding an electron at higher- 
energy near that of the peak. With decreased doping, 
the ratio of g(E) at negative and positive energies en- 
hances quite dramatically, e.g. from x=0.125 to 0.056, 
g(—A)/g(A) at the corresponding energy of the peak A 
(in units of t) increases from 1.96 to 2.73^. Similar be- 
haviors are found for the case with vanishing (t',t")/t 
(not shown). In contrast to this, for the d-BCS (inset in 
FigEJ) case in the same finite lattices there is almost no 
change of the ratio within the gap. The numerical results 
thus tells us features due to strong correlation which are 
not fully explored yet in the tunneling experiment, i.e. 
the particle-hole asymmetry of average conductance ex- 
ists even within the gap region and gets enhanced with 
underdoping. 

FigEI also reveals correlation between heights of the 
spectral weight peak and the gap size (or the width be- 
tween peaks) as doping level is varied 2 ^. Within the dop- 
ing level shown in Fig0J the peak height scales with the 
pairing amplitude but apparently anti-correlates with the 
gap size. This is in clear contrast to the BCS case in 
which the peak height, proportional to the SC coher- 
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FIG. 4: The doping dependence of the peak height for g(E), 
pairing amplitude and width between peaks in 20 x 20 lattice 
for d-RVB case. Data are extracted from that shown in Fig|^] 
and those obtained in the same way but not shown there. The 
values for the width between peaks are referred to the vertical 
axis on the right-hand side (indicated by the arrow) . 



ence, scales with the width between peaks or gap size as 
more holes doped into the system. Our result agrees 
qualitatively with what has been extracted from STS 
experiments^. 

To conclude, in order to provide a better understand- 
ing of the results measured by the tunneling experiments 
without the complication of mixing disorder and strong 
correlation, here we studied the SW's for adding and re- 



moving an electron for a uniform <i-RVB SC state without 
disorder. We derive analytically and examine numer- 
ically the relation between pairing amplitude and SW 
products. Performing particle-hole transformation, we 
obtain also exact dependence of SW for removing an 
electron with doping in the ED systems, which could be 
tested by photoemission spectroscopy. While the strong 
correlation effect is less noticeable by looking at the pair- 
ing amplitude, we found that the SW for removing an 
electron deviates clearly from results obtained by RMFT 
in the low doping regime. More specifically, at this dop- 
ing level the conductance-related quantity of the uni- 
form c?-RVB state on finite lattices computed exactly is 
particle-hole asymmetric below the gap energy and, con- 
sistent qualitatively with what is seen in recent tunneling 
experiments, the extracted gap value (from the excitation 
energy) or, equivalently, the width between SW peaks 
anti-correlates with the peak heights. 
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